
clear all
set more off

log using 01_analysis.log
use 01_main.dta


**============Table 1====================================
sum econ misconduct lifebehavior behavior

**============Figure 1====================================
tab chuliwithin if type==1
tab chuliwithin if type==2
tab chuliwithin if type==3
tab chuliwithin if type==4


tab discipline if type==1
tab discipline if type==2
tab discipline if type==3
tab discipline if type==4



**============Table 2====================================
sum tt_repost_ln chuliwithin discipline type rank_code position_code
sum econ misconduct lifebehavior behavior

tab rank_code,mi
tab position_code,mi
tab chuliwithin
tab discipline,mi

//Note: the detailed information on rainfall variable will be shown in 02_analysis.do

**============Table 3====================================
probit chuliwithin tt_repost_ln i.type,vce(robust)

probit chuliwithin tt_repost_ln i.type i.position_code i.rank_code,vce(robust)

oprobit discipline tt_repost_ln i.type,vce(robust)

oprobit discipline tt_repost_ln i.type i.position_code i.rank_code,vce(robust)




**============Table 4====================================

oprobit discipline tt_repost_ln jine_ln if type==1,vce(robust)	  
oprobit discipline tt_repost_ln jine_ln i.position_code i.rank_code if type==1,vce(robust)
	  

char   xingchouwen[omit] 2 
xi:oprobit discipline tt_repost_ln i.xingchouwen  if lifebehavior==1,vce(robust) 
xi:oprobit discipline tt_repost_ln i.xingchouwen i.position_code i.rank_code if lifebehavior==1,vce(robust)


**============Figure 2 [panel A]===========================

probit chuliwithin tt_repost_ln i.type i.position_code i.rank_code if rank_code<=5

oprobit discipline tt_repost_ln i.type i.position_code i.rank_code if rank_code<=5



**============Figure 2 [panel B]===========================


probit chuliwithin tt_repost_ln i.type i.position_code i.rank_code if  position_code==5 | position_code==6 | position_code==1 

oprobit discipline tt_repost_ln i.type i.position_code i.rank_code if  position_code==5 | position_code==6| position_code==1 



**============Table 5===========================

tobit discipline tt_repost_ln i.type i.position_code i.rank_code, ll(0) vce(robust) 

drop if discipline==0
oprobit discipline tt_repost_ln i.type i.position_code i.rank_code,vce(robust)


log close
